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Abstract. This paper presents the analysis of Doppler p-mode observations of the G2V star a Cen A obtained 
with the spectrograph CORALIE in May 2001. Thirteen nights of observations have made it possible to collect 
1850 radial velocity measurements with a standard deviation of about 1.5 ms~ . Twenty-eight oscillation modes 
have been identified in the power spectrum between 1.8 and 2.9 mHz with amplitudes in the range 12 to 44 
cms -1 . The average large and small spacing are respectively equal to 105.5 and 5.6 /^Hz. A comparison with 
stellar models of a Cen A is presented. 
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1. Introduction 

The lack of observational constraints leads to serious un- 
certainties in the modeling of stellar interiors and stellar 
evolution. The measurement and characterization of os- 
cillation modes in solar-like stars is an ideal tool to test 
models of stellar inner structure and theories of stellar 
evolution. 

The discovery of propagating sound waves, also called 
p- modes, in the Sun by Leigh ton et al. ( 1962] ) and their 
interpretation by Ulrich (1970), has opened a new area in 
stellar physics. Frequency and amplitude of each oscilla- 
tion mode depend on the physical conditions prevailing in 
the layers crossed by the waves and provide a powerful 
seismological tool. Helioseismology led to major revisions 
in the "standard model" of the Sun and provided for in- 
stance measures of the inner rotation of the Sun, the size 
of the convective zone and the structure of the external 
layers. 

Solar-like oscillation modes generate periodic motions 
of the stellar surface with periods in the range 3-30 min- 
utes but with extremely small amplitudes. Essentially two 
methods exist to detect such a motion: photometry and 
Doppler spectroscopy. In photometry, the oscillation am- 
plitudes of solar-like stars are in the range 2-30 ppm 
while they are in the range 10 - 150 cms -1 in radial 
velocity measurements. Photometric measurements made 
from the ground are strongly limited by scintillation noise. 
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To reach the needed accuracy requires observations made 
from space. In contrast, Doppler ground-based measure- 
ments have recently shown their capability to detect os- 
cillation modes in solar-like stars. 

The first good evidence of excess po wer du e to mode 
oscillations was obtained by Martic etal. ( 1999 ) on the F5 
subgiant Procyon. Bedding et al. ( 2001 ) obtained a quite 
similar excess power for the G2 subgiant (3 Hyi. These two 
results have independently been confirmed by our group 
based on observation s mad e with t he CORALIE spectro- 
graph (Carrier et al. |2001a| , |2001b[ ). 

A primary target for the search for p-mode oscillations 
is the solar twin a Cen A (HR5459). Several groups have 
already made thorough attempts to detect the signature 
of p-mode oscillations on this star. Two groups claimed 
mode detections with amplitudes 3.2- 6.4 greater than so- 
lar (Gelly et al. |1986[ Pottasch et al. |1992| ). However, this 
was refuted by three other groups who obtained upper 
limits of m ode a mplitudes of 1.4-3 times solar (Brown 
k Gi lliland |l990| ; Edmonds & Cram |l995t K jeldsen et al. 
1999 ). More recently Schou & Buzasi ( ^001 ) made photo- 
metric observations of a Cen A with the WIRE spacecraft 
and reported a possible detection. The first unambiguous 
observation of p- modes in this star was recently made with 
the spectrograph CORALIE mounted on the 1.2-m Swiss 
telescope at La Silla Observatory (Bouchy & Carrier 2001 
(hereafter referred as BC01); Carrier et al. 2001c ). We 
present here a revised and extended analysis of the acous- 
tic spectrum of this star and compare it with theoretical 
models. 
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Fig. 1. Radial velocity measurements of a Cen A. A 3 or- 
der polynomial fit subtraction was applied for each night. 
The dispersion reaches 1.53 ms' 1 . 

Table 1. Distribution and dispersion of Doppler measure- 
ments 



Date 


Nb spectra 


Nb hours 


a ( ms 1 ) 


2001/04/29 


220 


10.91 


1.58 


2001/04/30 


259 


11.36 


1.49 


2001/05/01 


268 


11.84 


1.39 


2001/05/02 


246 


11.71 


1.40 


2001/05/03 


267 


11.77 


1.53 


2001/05/04 


35 


1.43 


1.44 


2001/05/05 








2001/05/06 








2001/05/07 


206 


9.32 


1.84 


2001/05/08 


116 


4.77 


1.39 


2001/05/09 


96 


4.06 


1.61 


2001/05/10 


81 


3.30 


1.34 


2001/05/11 


56 


2.25 


1.88 



2. Observations and data reduction 

The observations of a Cen A were carried out with the 
CORALIE fiber-fed echelle spectrograph mounted on the 
1.2-m Swiss telescope at the ESO La Silla Observatory. 
A description of the spectrograph and the data reduction 
process is presented in BC01. 

a Cen A was observed over 13 nights in May 2001. A 
journal of these observations is given in Table 1. A total 
of 1850 optical spectra was collected with a typical signal- 
to-noise ratio in the range of 300-420 at 550 nm. Radial 
velocities were computed for each night relative to the 
highest-signal-to-noise-ratio optical spectrum obtained in 
the middle of the night (when the target had the highest 
elevation). The obtained velocities are shown in Fig. 1. 
The dispersion of these measurements reaches 1.53 ms -1 
and the individual value for each night is listed in Table 1. 

3. Acoustic spectrum analysis 




Frequency [mHz] 



Fig. 2. Power spectrum of the radial velocity measure- 
ments of a Cen A. 



3.1. Power spectrum 

In order to compute the power spectrum of the veloc- 
ity time series of Fig. 1 , we used the Lomb -Scargle mod- 
ified algorithm (Lomb |l976 , Scargle 1982|) for unevenly 
spaced data. The resulting LS periodogram, shown in 
Fig. 2, exhibits a series of peaks between 1.8 and 2.9 
mHz modulated by a broad envelope, which is the typ- 
ical signature of solar- like oscillations. This signature also 
appears in the power spectrum of each individual night. 
Toward the lowest frequencies (v < 0.6 mHz), the power 
rises and scales inversely with frequency squared as ex- 
pected for instrumental instabilities. The mean white 
noise level a ps , computed in the range 0.6-1.5 mHz, reaches 
2.39 x 10 -3 m 2 s -2 . Considering that this noise is gaus- 
sian, the mean noise level in the amplitude spectrum is 
4.3 cms -1 . With 1850 measurements, the velocity accu- 
racy corresponds thus to 1.05 ms -1 . 



3.2. Search for a comb-like pattern 

In solar-like stars, p-mode oscillations are expected to 
produce a characteristic comb-like structure in the power 
spectrum with mode frequencies v n ^ reasonably well ap- 
proximated by the simplified asymptotic relation (Tassoul 



1980) 



v n ,i « Av (n + - + e) - 1(1 + l)5v m /§ 

with Apq = < v n _\ — Vn-ij > and 
8vo2 = < v n ,o - Vn-1.2 > ■ 



(1) 



The two quantum numbers n and I correspond to the 
radial order and the angular degree of the modes, respec- 
tively. A^o and 5vo2 are the large and small spacing re- 
spectively. For stars of which the disk is not resolved, 
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Fig. 3. Autocorrelation of the power spectrum of a Cen A. 



only the lowest-degree modes (I < 3) can be detected. 
In case of stellar rotation, high-frequency p-modes need 
to be characterized by a third quantum number m called 
the azimuthal order: 



~ V n ,lfi + m/Prot 



(2) 



with —l<m<l and P ro t the averaged stellar rotational 
period. 

One technique, commonly used to search for period- 
icity in the power spectrum, is to compute its autocorre- 
lation. The result is shown in Fig. 3. The first identified 
peak at 5.7 /zHz corresponds to correlation between modes 
1=2 and 1=0, hence the small spacing <5^o2- The second 
peak at 11.5 /iHz corresponds to the daily alias. Peak 3 at 
50 /iHz is associated with correlation between modes 2=0 
and 1=1 and/or 1=1 and 1=2. The fourth peak at 55.8 /iHz 
is associated with correlation between modes 1=2 and 1=1 
and/or 1=1 and 1=0. The separation between peaks 3 and 
4 corresponds to the small spacing 6 vq%- Peaks 5 and 7, 
at 99.9 and 111.4 /xHz respectively, are associated with 
correlations between modes 1=0 and 1=2. Finally peak 6 
at 105.5 /iHz corresponds to the large spacing A^o. Other 
peaks correspond to correlations between modes and daily 
aliases. 



3.3. Echelle diagram 

In order to identify the degree I of each mode individually, 
the power spectrum has been cut into slices of 105.5 /iHz 
which are displayed above one another . Suc h a process 
was proposed and used by Grec et al. ( 1983 ) to identify 
the degree I of solar modes. The result, also called echelle 
diagram, is presented in Fig. 4. 

Fig. 5 corresponds to the merged echelle diagram, all 
the slices of 105.5 /iHz were summed up. 
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Fig. 4. Echelle diagram of a Cen A oscillations. Modes 
= 0, 1 and 2 and their aliases are identified. 
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Fig. 5. Sum of the echelle diagram of a Cen A oscillations. 
Modes I = 0, 1 and 2 and their aliases arc identified. 



Both Figs. 4 and 5 make it possible to identify modes 
of degree Z=0, 1 and 2 and their side-lobes due to the daily 
aliases at ± 11.57 /iHz. These two figures indicate that the 
modes follow an asymptotic relation reasonably well. 



3.4. Gap filling process 

In order to reduce the effect of single-site observation, the 
repetitive music process proposed by Fossat et al. (1999) 
was used. This method is based on the fact that the au- 
tocorrelation function of the full disk hclioscismological 
signal, after dropping quickly to zero in 20 or 30 minutes, 
shows secondary quasi periodic bumps, due to the quasi- 
periodicity of the peak distribution in the Fourier spec- 
trum. In the helioseismological signal, the first of these 
bumps appears at about 4.1 hours and corresponds to a 
A^o/2 = 67.5 /iHz periodicity of the modes. An obvious 
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Table 2. Mode frequencies (in /iHz) of a Cen A. 
Uncertainty on frequency is estimated to 0.46 /iHz. 



Fig. 6. Observational window response with (black line) 
and without (grey line) the gap-filling process. 



consequence is that simply replacing a gap by the signal 
collected 4.1 hours earlier or 4.1 hours later provides a 
very efficient gap filling method. It should be noted that 
such a process modulates the background noise with the 
67.5 /iHz periodicity, the fine adjustment of the method 
consists of placing the maxima of this modulation on the 
p-mode frequencies, so that any loss of information is lo- 
cated in the noise. This gap-filling process is not suitable 
if the modes do not follow an asymptotic relation. As pre- 
caution the frequency range should be cut in several parts 
and the fine tuning of the temporal periodicity adjusted 
in each of them. 

In the case of a Cen A, the periodicity of the modes is 
equal to Ai/q/2 = 52.8 /iHz which correspond to a period- 
icity in time of the signal of about 5.26 hours. Considering 
that a Cen A was observed during the longest nights for 
a little less than 12 hours, this gap-filling process makes 
it possible to fill the duty cycle up to 90 %. The power 
spectrum of a pure sinusoid using the observational win- 
dow of a Cen A with and without gap filling is presented 
in Fig. 6. Note that aliases at ± 1.35 /iHz are not reduced, 
due to the gap of the seventh and eighth nights during the 
run. 

The frequency range was cut into 3 parts (1.8 - 2.2 
mHz, 2.15 - 2.55 mHz and 2.5 - 2.9 mHz) and the tempo- 
ral periodicity was adjusted in each of them to place the 
maxima of the modulation on modes 1=1 which are easily 
identifiable. For the 3 parts of the power spectrum, the 
temporal periodicity used was respectively 5.299, 5.294 
and 5.289 hours. The result is presented in Fig. 7. 

This gap-filling process does not cause a drastic change 
in the power spectrum below 2.5 mHz. Above 2.5 mHz, 
the process helps in the identification of modes, especially 
modes 1=0 and 1=2 by removing or reducing aliases. Note 
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2679.8 


2733.2 
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2786.2 


2837.6 


2887.7 


A v 




105.6 ± 0.5 


105.6 ± 0.4 


105.0 ± 0.5 



that this process is based on the assumption that frequen- 
cies follow an asymptotic relation and can destroy any 
other frequency and/or enhance their side-lobes. 

3.5. Frequencies of p-modes 

The strongest peaks of the power spectrum computed with 
the gap filling process were identified and are listed in 
Table 2. These peaks were selected such as to have a signal- 
to-noise ratio in the amplitude spectrum greater than 3 
(0.017 m 2 s~ 2 in the power spectrum). 

The frequency resolution of our time series is 0.93 /iHz. 
If we suppose that no mode is resolved, the estimated 
uncertainty on the frequency determination is 0.46 /iHz. 
The values of n and I are deduced from the asymptotic 
relation (see Eq. (1)) assuming that the parameter e is 
near the solar value (e Q ~ 1.5). Our first identification in 
BC01 based on the 5 first nights did not make a distinction 
between the modes 1=0 and 1=2. The last row in Table 2 
gives the average large spacing for each angular degree. 
All the identified modes in Table 2 are indicated in Fig. 7. 

The average large and small spacing and the parameter 
e are deduced from a least-squares fit of Eq. (1) with the 
frequencies of Table 2: 

A^ = 105.5 ± 0.1 (iHz, 

5vq2 = 5.6 ± 0.7 iiHz, 

e = 1.40 ± 0.02 . 

Fig. 8 represents the echelle diagram of peaks listed in 
Table 2. 



3.6. Amplitudes of p-modes 

In order to estimate the average amplitude of the modes, 
we followed the method of Kjeldsen & Bedding (1995). 
A simulated time series was generated with artificial sig- 
nal plus noise using our window function. The generated 
signal corresponds to modes centered on 2.36 mHz, sep- 
arated by Afg/2 and convolved by a gaussian envelope. 
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Fig. 7. Power spectrum of a Cen A with (black line) and without (grey line) the gap filling process. The identified 
frequencies are indicated with dotted lines. 
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Fig. 8. Echelle diagram of identified modes in the power 
spectrum. The modes 1=2 (a), 1=0 (□) and 1=1 (■) are 

represented from left to right. 



g 30 



— 
- 

£ 20 
o. 
E 
< 




2,2 2,4 2,6 

Frequency [mHz] 



Fig. 9. Doppler amplitude of identified modes of degree 
1=0 (□), 1=1 (■) and 1=2 (a). 



The generated noise is gaussian with a standard deviation 
of 1 ms -1 . Several simulations show that the envelope 
amplitude is in the range 29-33 cms" 1 . 

Table 3 and Fig. 9 present the estimated Doppler am- 
plitude of each identified mode with the assumption that 
none of them is resolved. This amplitude was determined 
as the height of the peak in the power spectrum after 
quadratic subtraction of the mean noise level. Fig. 9, like 
Fig. 5, shows that the modes 1=0, 1 and 2 have a similar 
average amplitude. 

In order to compare the observed amplitude with the 
expected one, we computed the amplitude ratios between 
mode /=0 and the other modes (/, m). These amplitude ra- 
tios (see Table 4) were determined following Gouttebroze 
& Toutain (1994) assuming that all the modes have the 
same intrinsic amplitude. The effect of limb darkening was 
neglected and the inclination of the rotational axis was as- 



Table 3. Doppler amplitude of identified modes deter- 
mined as the height of the peak in the power spectrum 
after quadratic subtraction of the mean noise level. 
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sumcd equal to the inclination of the orbital axis of the 
binary system (79°) determined by Pourbaix et al. (1999). 

The observed amplitude of modes 1=2 is higher than 
expected; this can be due to interferences between modes 
of different m. The modes corresponding to 1=3 are ex- 
pected to have an amplitude of only 20 % of the 1=1 mode, 
hence to be well below the threshold of identification. 
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Table 5. Observational constraints and global character- 
istics of the models. 



-2-10 1 
Frequency [/iHz] 

Fig. 10. Superposition of identified modes 1=0 (□), 1=1 
(■) and 1=2 (a) in the power spectrum with frequency 
shifted to zero. 



3.7. Linewidth of p-modes and rotational splitting 

The expected linewidth for a 1.1 M© star in the fre- 
quency range 1.8 to 2.9 mHz is less than 2 /iHz (Houdek 
1999). Our frequency resolution of 0.93 /iHz is then too 



low to attempt a clear determination of the linewidth of 
the identified modes. Below 2.2 mHz, modes seem un- 
resolved and present a linewidth equal to the frequency 
resolution, hence a damping time longer than 25 days. 
Above 2.2 mHz, the structure of modes seems more com- 
plex. This could be eventually caused by the fact that 
their lincwidths begin to be greater than the frequency 
resolution. 

For a Cen A, Hallam et al. ( |l99l[ ) estimated a ro- 
tational p eriod of 28.8 ± 2.5 days. More recently Saar 
& Osten ( 1997 ) measured a rotational velocity of 2.7 ± 
0.7 kms" 1 . With an estimated radius of 1.2 R Q , the pe- 
riod of rotation is 22 days. Assuming a uniform rotation, 
the corresponding splitting of the modes is expected to be 
about 0.5 /iHz. Assuming an inclination of the rotational 
axis of 79°, a separation of about 1 /iHz should be visible 
between modes (i=l,m=-l) and (/=l,m=l), and modes 
(Z=2,m=0) and (Z=2,m=±2). 

Figure 10 shows the superposition in the power spec- 
trum of modes 1=0, 1=1 and 1=2 with the frequency of 
each mode shifted to zero. Modes /=0 present side lobes 
at ± 1.35 /zHz, probably introduced by the aliases due to 
the gap of the seventh and eighth nights during the run. 
Modes 1=2 present side lobes at ± 1.25 /iHz which could 
be introduced by these same aliases or/and a rotational 
splitting. In the latter case it should correspond to a pe- 
riod of rotation of about 19 days, lower than the expected 



Models 


Morel et al. 


Guenther & Demarque 


M/M Q 


1.16 ±0.031 


1.1015 ±0.008 


T off [K] 


5790 ± 30 


5770 ± 50 


logs 


4.32 ± 0.05 


4.28 ± 0.02 




0.20 ±0.02 


0.22 ±0.02 


L/Lq 


1.534 ±0.103 


1.572 ±0.135 


i[Myr] 


2710 


5640 


Yi 


0.284 


0.300 


(f)i 


0.0443 


0.0480 


a 


1.53 


1.86 



one. Modes 1=1, which show a side lobe only on the left, do 
not allow a conclusion. The CLEAN algorithm (Roberts et 
al.~ 



1987 ) was used in order to remove the effect of the ob- 
servational window response. We noticed that this process 
docs not eliminate the whole signal at modes 1=2, as is the 
case for 1=0 and 1=1. The pattern shown in Fig. 10 leads 
one to think that an error of ± 1.3 /iHz could have been 
introduced at some identified mode frequencies and could 
explain the dispersion of the mode frequencies around the 
asymptotic relation shown in Fig. 8. 

Additional data will be needed to determine both the 
rotational splitting and the damping time of a Cen A p- 
modes. 



4. Comparison with models 

Seismological parameters deduced from our observations 
are in full agreement with the expected values scaling from 



the Sun (Kjeldsen & Bedding, 1995) giving the frequency 
of the greatest mode ^ max = 2.3 mHz, the large spac- 
ing Ai>q = 105.8 /iHz and the oscillation peak amplitude 
A>sc = 31.1 cms" 1 . 

The variations with frequency of the large and small 
spacings are compared in Figs. 11 and 12 with the val- 
ues deduced from models recently developed by Morel et 
al. ( [2000D and Guenther & Demarque ( |2000| ). The observa- 
tional constraints for the model calibrations and the prop- 
erties of the calibrated models are given in Table 5. The 
two sets of models differ essentially by the estimated mass 
and the age of the star. The large spacing, mainly related 
to stellar density, varies proportional to M 1 / 2 R~ z / 2 . The 
small spacing, mainly related to the structure of the core, 
decreases with stellar age and mass. 

Our observations are well bracketed by the two sets of 
models for both the large and small spacing. This clearly 
suggests a model with intermediate properties, i.e. a mass 
between 1.10 and 1.16 Mq. New mo dels of a Cen A are 
being developed by Thevenin et al. ( [2002 ), with different 
properties in order to better interpret the observed oscil- 
lations and constraint the physical parameters of this star. 
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Fig. 11. Variations of the large spacing between modes 
with consecutive radial order A^ n ,; = v n ,i — v n -\.i for p- 
modes of degree ^=0 (□), 1=1 (■) and 1=2 (A) versus 
frequency. Upper lines correspond to the model of Morel 
et al. and lower lines to that of Guenther & Demarque. 
Solid, dashed and dotted lines correspond respectively to 
modes ^=0, 1 and 2. 



5. Conclusions and prospects 

Our observations of a Cen A yield a clear detection of p- 
modc oscillations. Several identifiable modes appear in the 
power spectrum between 1.8 and 2.9 mHz with an average 
large spacing of 105.5 /iHz, an average small spacing of 
5.6 /iHz and an envelope amplitude of about 31 cms -1 . 
These characteristics, in full agreement with the expected 
values, make it possible to put constraints on the physical 
parameters of this star. 

Additional data with a higher signal-to-noise and a 
higher frequency resolution will make it possible to deter- 
mine properly the suspected rotational splitting and the 
damping time of a Cen A p-modes. 

This result, obtained with a small telescope, demon- 
strates the power of Doppler ground-based asteroseismol- 



ogy. Future spectrographs like HARPS (Queloz et al. 2001 



Bouchy 2001) are expected to conduct asteroseismological 
study on a large sample of solar- like stars, and to enlarge 
significantly our understanding of stellar physics. 
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Fig. 12. Variations of the small frequency spacing 5vq2 = 
v n Q — v n -\,2 versus frequency. The upper line corresponds 
to the model of Morel et al. and the lower line to that of 
Guenther & Demarque. 
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